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Abstract 

A formulation of non-relativistic quantum mechanics in terms of Newtonian 
particles is presented in the shape of a set of three postulates. In this new 
theory, quantum systems are described by ensembles of signed particles which 
behave as field-less classical objects which carry a negative or positive sign 
and interact with an external potential by means of creation and annihila¬ 
tion events only. This approach is shown to be a generalization of the signed 
particle Wigner Monte Carlo method which reconstructs the time-dependent 
Wigner quasi-distribution function of a system and, therefore, the correspond¬ 
ing Schrodinger time-dependent wave-function. Its classical limit is discussed 
and a physical interpretation, based on experimental evidences coming from 
quantum tomography, is suggested. Moreover, in order to show the advantages 
brought by this novel formulation, a straightforward extension to relativistic ef¬ 
fects is discussed. To conclude, quantum tunnelling numerical experiments are 
performed to show the validity of the suggested approach. 

Keywords: Mathematical formulation of quantum mechanics, Newtonian 
particles, Wigner formalism, Schrodinger equation. Relativistic quantum 
mechanics. Quantum tunnelling 


1. Introduction 

Nowadays, many different mathematical formulations of quantum mechan¬ 
ics exist, among which the ones suggested by E. Schrodinger [T], E. Wigner [2], 
R. Feynman [3], L.V. Keldysh [3], K. Husimi [^, D. Bohm m, m are the most 
popular ones. While, at a first glance they seem to be drastically different the¬ 
ories, it can be shown that they are all mathematically equivalent. As a matter 
of fact, they offer the same set of predictions and can be seen as complementary 
points of view. The situation is not any different than what one observes in 
classical mechanics where different formalisms, such as Newtonian, Lagrangian, 
Hamiltonian, etc., coexist and shed different light on different mechanical as¬ 
pects. 
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In this paper, we introduce a new formulation of quantum mechanics in 
terms of signed classical field-less particles. This suggested theory is based on a 
generalization and a novel physical interpretation of the mathematical Wigner 
Monte Carlo method [5] which is able to simulate the time-dependent single- 
and manj-body Wigner equation 0, m in a quite intuitive fashion, which 
experimentalists are familiar with. Indeed it describes quantum objects in terms 
of classical particles only. 

One should notice that the signed particle formulation is equivalent to the 
usual formulations. As such, no fundamental new results are introduced. The 
predictions made are the same as the ones made in the more standard theory. 
However, ’’there is a pleasure in recognizing old things from a new point of 
view” [5] and the author hopes it can offer a new perspective on the puzzling 
quantum nature of the microscopic world. The new theory is based on clas¬ 
sical particles which have a position and momentum simultaneously although 
the Heisenberg principle of uncertainty is still respected in the formulation and 
embedded in terms of unknown initial conditions. In particular, we will show 
that the sign of a particle cannot be evaluated experimentally and no physical 
measurement can be depicted to find differences with other formalisms. Nev¬ 
ertheless, it offers some noticeable advantages. On the one hand it is a very 
intuitive approach which provides a new way to describe Nature at a quantum 
level. On the other hand, it is a computationally attractive formulation being 
based on independently evolving particles, allowing deep levels of parallelization 
in the time-dependent simulation of quantum single- and many-body systems. 
Finally, it allows the inclusion of physical effects which are difficult to treat in 
other formulations of quantum mechanics. 

In this paper, we first introduce the three postulates which completely de¬ 
fine the new mathematical formulation of quantum mechanics in terms of signed 
particles. We then proceed with showing that these three postulates are enough 
to replicate the results of more conventional quantum theories. In particular we 
show that our suggested approach is a natural generalization of the Wigner 
Monte Carlo method which reconstructs the time-dependent Wigner quasi¬ 
distribution function and, thus, the corresponding Schrodinger wave-function. 
Afterwards, its classical limit is considered in several details. In order to show 
the applicability of the proposed theory, we numerically simulate several exper¬ 
iments involving quantum tunnelling. To conclude, we extend the theory to 
include general relativity in the formalism by generalizing the second postulate 
to the case of particles travelling along space-time geodesics, and we show that 
typical relativistic effects are observable for quantum wave-packets. 

2. A signed particle approach 

We now introduce the signed particle formulation of quantum mechanics by 
enunciating, and discussing, the three postulates which completely define the 
new theory. 
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Postulate I. Physical systems can be described by means of (virtual) New¬ 
tonian particles, i.e. provided with a position x and a momentum p simultane¬ 
ously, which carry a sign which can be positive or negative. 

Postulate II. A signed particle, evolving in a potential V = V (x), behaves 
as a Geld-less classical point-particle which, during the time interval dt, creates 
a new pair of signed particles with a probability 7 (x(t)) dt where 



( 1 ) 


and V(( (x; p) is the positive part of the quantity 



known as the Wigner kernel (in a d-dimensional space) |3- moment 

of creation, the parent particle has sign s, position x and momentum p, the 
new particles are both located in x, have signs +s and —s, and momenta p + p^ 
and p — p' respectively, with p' chosen randomly according to the (normalized) 


probability 


Postulate III. Two particles with opposite sign and same phase-space co¬ 
ordinates (x, p) annihilate. 

We will show that this set of three postulates is sufficient to reconstruct 
the time-dependent quasi-distribution function of a system and, consequentely, 
its wave-function. In other words, the time-dependent evolution of a quantum 
system can be completely expressed in terms of creation and annihilation of 
signed particles only. Some comments on the interpretation and application of 
the postulates follow below. 

Momentum integral. It is important to note that the definition of the func¬ 
tion 7 = 7 (x) is introduced in a continuum phase-space which, in turn, intro¬ 
duces a non-Riemann integral Q. This integral is, in the spirit, very similar to 
a path integral and the author proposes to name it a momentum integral. In 
fact, while a path integral is spanning the space of trajectories, the integral 0 
similarly spans the space of momenta. Of course, this new definition is rather 
elementary and necessitates more mathematical investigations as it may not 
always imply convergence. 

Particles with a sign. The concept of signed particles was suggested for the 
first time in m to explain the emergence of negative peaks in the Wigner quasi¬ 
distribution of particular quantum systems (and the reader should note that it 
was introduced, strictly speaking, as a mathematical tool). Consecutively, this 
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concept has been applied with success to the numerical simulation of a challeng¬ 
ing quantum benchmark test involving the evolution of 5—functions |12j . Thus, 
we borrow this notion to enunciate the first postulate. In particular, positive 
and negative particles can be interpreted as usual classical field-less particles 
moving on Newtonian trajectories, but it is impossible to experimentally know 
simultaneously their position and momentum due to the Heisenberg uncertainty 
principle, thus denying the chance of calculating their trajectory (and one can 
only start from a distribution of particles in practical calculations). The physi¬ 
cal interpretation of the negative sign is suggested below based on experimental 
evidences coming from the field of quantum tomography a set of 

techniques to reconstruct the Wigner function of a given experimental setting. 

Negatively signed particles. In practice, an experiment is prepared in spe¬ 
cific initial conditions and repeated a large amount of times (quantum systems 
have different outcomes every time), providing a projection of the correspond¬ 
ing quasi-distribution on the spatial coordinate axis and momentum axis, and 
by applying the inverse Radon transformation one reconstructs a higher dimen¬ 
sional function which is the Wigner function of the system (with positive and 
negative values). Now, classical objects are always localized in a precise point 
of the phase-space (represented by <5—functions), while for quantum objects the 
presence of the Heisenberg principle of uncertainty prevents such localization, 
forcing the description of the dynamics to an area of the phase-space bigger than 
AxAp = |. This is clearly exhibited by the appearance of negative values in 
the Wigner function and, therefore, one may infer that particles with a negative 
sign are essentially those which are experimentally unreachable El- 

Physical picture. Consequently, the picture that this novel theory offers is 
rather peculiar and different than any other mathematical formulation of quan¬ 
tum mechanics. Quantum systems are now described by means of ensembles 
of Newtonian field-less particles which now carry a sign and interact with an 
external potential by means of creation and annihilation events only. When a 
pair of particles is created, one is in an experimentally reachable state (positive 
sign), and the other in a non-reachable one (negative sign). This offers a phys¬ 
ical picture which is relatively easy to grasp and which allows the inclusion of 
quite complex effects in a natural way (as it will be shown below). 

Practical implementation of the formulation. In particular, given a signed 
particle at time t with sign s, mass m, and (phase-space) coordinates (x, p), 
indicated from now on by (s,m;x, p), we introduce the operator S which con¬ 
structs a new set of three signed particles 

S[{s,m;x,p)] = {(s,to;xi,pi) , (-ks,m;X 2 ,P 2 ), (-s,m; X3, pa)} 
in the following way: 

• At time t, one generates a random number r S [0,1] and computes the 
quantity St = 
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• At time t + St, the initial particle evolves as field-less and has new coordi¬ 
nates (xi,pi) = (^x-H p^, 

• A pair of new signed particles is created at time t + St, where the par¬ 
ticle with the sign s has coordinates (x 2 ,P 2 ) = ^x-|- ^^^t,p-|-p'^ and 

the particle —s has coordinates (xsjPs) = - 1 - ^^St,p — p'^ and the 

quantity p' is computed from the normalized probability 

The reason for generating a random number r is easily explained. The 
probability that a signed particle, in a small interval dt, generates a new pair of 
particles is 7 (x(t)) dt. Consequently, the probability that a particle which has 
generated a new signed pair at time to has not yet generated another new pair 
at time t is 

- /,* dd7(x(t')) 
e •’to ’^ 

or, in other words, the probability that no creation occurs during the interval 
[to,t]. Therefore, the probability P{t) that a particle will generate a new pair 
during the interval dt is given by 

P[t)dt = ^ {-x.{t)) e dt. 

Due to the complexity of evaluating the integral in the exponent, we thus utilize 
the numerical technique developed in the context of the Boltzmann Monte Carlo 
method for phonon scattering events m, which consists of generating a random 
number r in the way described above to find when a (creation) event happens. 

In practical calculations, one starts from an ensemble of particles distributed 
in the phase-space according to some specified initial quasi-distribution function, 
which may be obtained, for example, from a particular wave-function or density 
matrix describing the initial conditions of the experiment (but not restricted 
to these mathematical objects only). It is trivial to generalize the operator S 
to act on a set of N signed particles (with N an arbitrary natural number) by 
applying the above algorithm to every single particle in the given set. Further¬ 
more, if a time interval [t, t + At] is provided, the operator S can be applied 
repeatedly to every particle appearing in the creation process until the final time 
t -I- At is reached. The final set of particles is, therefore, constituted of every 
particle involved in the process described above. It is also important to men¬ 
tion that practical use of the third postulate makes sense only when the time 
and phase-space are discretized and two signed particles annihilate when they 
belong to the same cell and have opposite signs at the same time step. Thus, 
for computational purposes involving the postulates II and III, discretized time 
and phase-space have to be introduced (although this does not preclude the 
possibility of an analytical/theoretical treatment on a continuum phase-space 
in certain circumstances). 

For the sake of clarity, a schematic of the algorithm is presented below (the 
reader is also encouraged to download the working source code from [E]). 
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1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 
9. 
10 


start from an initial ensemble of signed particles 
compute the Wigner kernel 
compute the gamma function 
for time from 0 to final time 

for every particle of the ensemble 
compute free flight 

(eventually) create new pairs of signed particles 
end for 
annihilation 
end for 


In particular, row 1 is equivalent to assigning initial conditions to the system, 
row 2 is based on formula <§ - eventually restricted to a finite and discrete 
domain -, row 3 is obtained by means of Q - eventually restricted to a finite and 
discrete domain -, and, finally, the nested loops starting on row 5 are equivalent 
to the recursive application of the operator S. Moreover, the reader should 
note that row 9 suggests the application of an annihilation step. This is a very 
important point as, in certain situations, the number of signed particles may 
increase indefinitely if no annihilation scheme is applied (therefore corrupting 
the validity of the solution). When exactly to apply this technique during the 
simulation remains an open question which will be addressed elsewhere. 


On macroscopic variables. The signed particle formulation makes predictions 
by averaging a macroscopic variable ^ (x; p) over an initial set of particles, 
taking into account their sign, i.e. 


1 « , 
i=l 

where si, and Pi are the sign, the position and the momentum, respectively, 

of the i-th particle. Some well-known examples of macroscopic variable are 

2 

A (x; p) = 1 which corresponds to the probability density, A (x; p) = ^ which 
corresponds to the kinetic energy, etc. This is similar to what is done in classical 
statistical mechanics within the context of the Boltzmann Monte Carlo method 
m- Given a function A = >1 (x; p), formula ([^ represents the only meaningful 
way to make predictions in this theory and which can be compared to experi¬ 
mental measurements. Finally, one should note that the sum in ([^ runs over 
the whole set of N particles. Any other use involving subsets of signed particles 
is not allowed in the theory. 

3. Derivation of the postulates 

In this section, we show that the set of three postulates defined in the previ¬ 
ous section has its origins in the time-dependent single-body Wigner equation 
[2]. In fact, algebraic manipulations of this equation endorse a unique physical 
interpretation which, eventually, offers a novel physical picture on the dynamics 


6 


of quantum systems. In a sense, this is not any different from other formulations 
which come from previous formulations (e.g. the Wigner formalism is directly 
obtained from the time-dependent Schrodinger equation, although it allows a 
totally different perspective). Finally we show that, in this context, the signed 
particle Wigner Monte Carlo (MC) method (defined and thoroughly described 
in 0, i, nni) can be naturally interpreted as a numerical discretization of the 
new formulation (by introducing a semi-discrete phase-space), restricted to a 
finite domain. 

The time-dependent single-body Wigner equation [5] describes the evolution 
of a quasi-distribution function fw = /w(x; p; t) defined over a phase-space and 
reads: 

^^(x;p;t)-t — • Vx/vv = / dp7w(x; p-H p'; t)IV(x; p'; t), (4) 

at m J 

with Vw = lV(x;p;i) the Wigner kernel (§ m- It is possible to show that 
the prediction made by means of this equation are exactly the same as the 
ones provided by the Schrodinger equation (through the invertible Wigner-Weyl 
transform) [T^]. Therefore, solving the Wigner equation is equivalent to solve 
the Schrodinger equation, despite they provide two very different pictures. 

Through simple algebraic manipulations (by adding to both sides of the 
equation the quantity 7 (x(?/)) and reordering the terms), it is possible to ex¬ 
press the Wigner equation @ in a integral form, more precisely as a Fredholm 
equation of second kind: 

/H.(x;p;f)-e-/o^«^»‘^V.(x(0);p) = 

/*oo ^“i-OO p p± 

/ dt' / dp' / dx'/^y(x'; p'; t')r(x'; p; p')e~ Jt'— t')(5(x'— 
Jo J — oo J 


where 


r(x, p, p') = V+ (x, p - p') - V+ (x, -(p - p')) + 7 (x)d(p - p'), 

with fi (x; p) formally representing the initial conditions of the system at time 
t = 0. The reader should note that no discretization of the phase-space is 
introduced in the process. 

Now, equation ([^ suggests that it is possible to express the expectation 
value (A) of a macroscopic variable A = A (x; p) as a Liouville-Neumann series. 
For instance, the zeroth order term of the series reads: 

poo p p-\-CC> nt' 

(^)o(t) = / dt' dx dp'/j(x;p')e" Jo '^(^(*'^^''^A(x(t'),p')d(t'- r). 

By simple mathematical considerations, it is possible to interpret part of the 
integrand as a product of conditional probabilities. Thus, if fi = /i(x;p) is 
normalized to unity, one generates random points (x; p') in the (continuum) 


(5) 

x(t'))6»£,(x'), 
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phase-space at the initial time t = 0. This process initializes the particle trajec¬ 
tories X = x(t) and the exponent gives the probability for the particle to remain 
over the trajectory provided that the out-of-trajectory event rate is provided 
by the function 7 = 7(x). This probability filters out these particles, so that 
the randomly generated out-of-trajectory time is less than r. If the particle re¬ 
mains in the same trajectory till time r, it has a contribution to (^)o(t) equal to 
/i(x, p')A(x(t), p'), otherwise it does not contribute at all. Therefore, (A)o(r) 
is estimated by the mean value obtained from the N initialized particles. 

In the same way, the first three terms of the Liouville-Neumann series show 
how to continue with higher order terms in a way similar to [5] (but over a con¬ 
tinuum phase-space). Indeed, one observes that the expansion of (A) actually 
branches, and the total value is given by the sum of all branches. Considering 
the presence of the term T = r(x, p, p') in the series, one can equivalently talk 
in terms of three appearing particles. The initial parent particle survives and 
a couple of new signed particles, one positive and one negative, are generated 
with the first two probabilities. This introduces the first postulate of our sug¬ 
gested formulation (signed particles) and part of the second postulate (particles 
generation). 

In the creation process, one momentum state p — p' = q is generated with 
probability: 

7 (x) 

and, with the same probability, we generate another value, say q', for the second 
momentum state p' — p = qb Thus, after any free flight the initial particle 
creates two new particles with opposite signs and momentum offset (around 
the initial momentum) equal to -l-q and — q with q = p — pb This completes 
the formulation of the second postulate. The third and last postulate (particles 
annihilation) is a simple consequence of the way one computes macroscopic 
variables, as depicted in formula (|^. 

This concludes the derivation of the set of postulates presented in the pre¬ 
vious section. In this context, it now becomes clear that, if one restricts the 
spatial domain over a finite length and introduce a discretization of the momen¬ 
tum space, the Wigner MC method is re-obtained m,m- Therefore, the signed 
particle formulation of quantum mechanics represents a generalization of this 
method which can be considered a valid new physical picture of the quantum 
world. 


4. The classical limit 

In this section we explore the classical limit of our suggested formulation 
for > 0. To this aim, we distinguish between two cases: dynamics in free 
space and dynamics in the presence of an external potential. This distinction 
allows a natural and simple way to investigate how the signed particle theory 
can re-obtain the Newtonian theory. Although many assumptions are made in 



this section, we obtain an important point: signed particle creation is prohibited 
in the classical limit. 


Free space. In the case of free space (absence of any external potential), and 
in the classical limit 7i —>■ 0, the Heinsenberg principle does not represent a re¬ 
striction any longer. Therefore, it is possible to start from an initial distribution 
function for the signed particles which has the shape of a Dirac <5—function (in 
other words, one particle with a given position and momentum). Furthermore, 
in the classical limit, particles have always a positive sign as it is well-known 
that classical distribution functions are positive definite. In this context, thus, a 
signed particle behaves as a classical Newtonian (field-less) particle as expected 
and the classical limit is trivially obtained. It is interesting to note that in the 
absence of any external potential the function 7 (x) is identically zero or, in 
other words, there is no creation of particles. 


External potential. We now discuss the case in the presence of a non- 
identically-zero external potential V = V(x) in the classical limit h —)■ 0. For 
the sake of simplicity, we restrict ourselves to the one-dimensional case, as the 
extension to the three-dimensional space is trivial. 

One starts by noting that given a potential, and having calculated the func¬ 
tion 7 ( 0 :) from it, the probability for a (positively signed) particle to stay in a 
given path x = x(t), during the time interval [ 0 ,t], is 


Now, one can show that the Wigner kernel ([^ reads 


^ + 00 


Vwix;p) = —2 / 

TT/i J — Q 


dt' 


_P±C_ p 

e — 


vL + -v(x- 

\ 2m / V 2m / 


under the change of integration variable x' = ^t'. By exploiting the fact 
that the Wigner kernel is a real function [2] and supposing that the potential 
V = V(x) can be expressed in terms of a McLaurin series, one easily obtains 
(by means of trivial algebraic manipulations) the following expression for the 
kernel 


Vw{x;p) 


2 ^ /■+“,, ,pH\V^'^\x) 

/t—0,et;en ^ 

rdt'cos Avix)F 

7-00 

p^t' 1 d'^V(x) f pt' 


nh^m 


2p 

nh^m J-ao 

t4 




v^t' 

/ dt'cosA —) 

7-00 nm 


2 dx^ 


2m 



n 


2 

-f 


where the summation runs over even numbers only and (x) is the n—th 
order derivative of the potential. Finally, by neglecting every term down to the 
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second order 0(h?), and by performing the analytical integration of the only 
term remaining, one has 


Vw{x;p) 


2V{x) . . 

nph 


2,1 -1 +00 
P^t 


— oo 
2 +/ 


2V(x) , 

— -— 

t'-^+oo TTpn Tim 


2 lim 

t'-^+c 

W^x) 


sm( —), 
Trpfi, Tim 


) 


(for a big enough time T) which substituted in 0 gives the following expression 
for the function 7 ( 0 ;) 


7 ( 0 :) = lim 
Ap-)-0+ 


+ 00 

E 

M——00 


w(x) . AM/T^pfr: 

■sin(-V--) 


TTTi{MAp) 


Tim 


with [.. .]^ the positive part of the quantity in brackets. Now, supposing that 
V{x) > 0 for any position x (the extension to a general potential is trivial), the 
condition for the term in the bracket to be positive reads: 


0 < M < 


Timir 

Ap2r’ 


which, in the classical limit ?i — 0 , becomes 


0 < M < 0, 

or in other words M — 0. Therefore (in the classical limit) the function 7 is 
identically zero 7 (x) = 0 , and creation of pairs of signed particles is prohibited. 

Finally, concerning the equations of motion for (positively) signed particles 
for ?i —0, one starts from the fact that the Wigner equation reduces to the 
Vlasov equation [5]. Therefore, the technique utilized in the previous section to 
obtain the postulates of our formulation must be applied to the Vlasov equation. 
By exploiting the main results of |17j (which is an application of the same tech¬ 
nique to the Boltzmann equation) one easily recover the Newtonian equations. 
This completes the study of the classical limit. 


5. Simulation of quantum tunnelling 

Tunnelling is a typical quantum effect which cannot be explained in terms of 
classical mechanics. In fact, it represents one of the foundational experimental 
evidence for the need of a quantum theory. In particular, it shows that material 
particles can tunnel through potential barriers, even when the initial particle 
energy is classically unsufficient. In this section, we show that such experimen¬ 
tal observation can be reproduced by our suggested formulation of quantum 
mechanics, thus proving the validity and the applicability of the approach. 
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Two benchmark tests are performed (similar to the one presented in [S]) 
consisting of an initially Gaussian wave-packet with energy equal to 0.08 eV 
directed towards a potential barrier equal to 0.10 eV (in the first case) and 
to 0.30 eV (in the second case). The simulation of a free wave-packet is also 
performed for comparison purposes and shown in Fig. In more details, in the 
first numerical experiment (0.10 eV) one expects that part of the wave packet 
to be reflected back and part of the packet to tunnel through the barrier (due 
to its energy slightly smaller than the barrier itself). In the second experiment, 
instead, a strong reflection is expected due to the high energetic value of the 
barrier (0.30 eV). 

Every simulation is fully obtained by evolving an initial set of signed particles 
as prescribed in the three postulates (introduced above) and augmented by the 
equations of motion for a signed particle: 

x(t) = x(t — At) + —--At 

m 

p{t) = p{t-At). (6) 

The results of these two experiments are reported in Figs. and showing 
the simulation results at times (from left to right, from top to bottom) 0 fs, 10 
fs, 20 fs, 30 fs, 40 fs, 50 fs, 60 fs and 70fs respectively. In particular, one notes 
in both experiments the appearance of negative values during the evolution, 
a typical sign of the presence of dominant quantum effects (in other words 
tunnelling). Furthermore, it becomes clear that a certain amount of signed 
particles tunnels through the barrier (represented schematically by a red line in 
the middle of the domain) although this amount varies according to the height 
of the barrier. We believe that the fact that quantum tunnelling effects can be 
reproduced by applying the operator S shows the validity and applicability of 
our suggested formulation of quantum mechanics. 

6. Relativistic formulation 

An extension to relativistic quantum mechanics can be delineated taking 
into account the classical nature of the particles involved in the new formula¬ 
tion of quantum mechanics. Indeed, when introducing the second postulate, 
no equation for the trajectory is specified. The enunciation is very general and 
expressed only in terms of classical Beld-less particles. Eventually, the equa¬ 
tions are specified for the first time in ([^ when necessary for implementation 
purposes, but this is not the only option available. 

As a matter of fact, it is possible to define the trajectories of particles in¬ 
volved in postulate II as geodesics, allowing the evolution of quantum objects in 
a curved space-time. Hence, given a metric tensor describing the geomet¬ 
ric structure of space-time, one can substitute eqns. by a relativistic orbit 
(geodesics): 
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where a, /3, /j, = 0... 3, s is the proper time, are the Christoffel symbols (a 
function of and = t, x^ = x, x^ = y, x^ = z (x, y, z spatial coordinates). 
In particular, if the metric tensor describes a flat space-time (Minkowski metric), 
one obtains special relativity. 

Figs. 1^ and show the results obtained from the simulation of an initially 
Gaussian wave-packet moving in a flat space-time with no external potential at 
a speed comparable to the speed of light (about |c). In particular. Fig. |^shows 
the time dilation computed (as a macroscopic average over all involved particles, 
see (§) in the non-relativistic (blue dashed curve) and relativistic (red contni- 
uous curve) cases respectively. Fig. instead, displays the (damped) induced 
dipole from which one obtains the energy spectrum of the system through a 
Fourier transform. It is clear from these pictures how typical relativistic effects, 
such as time dilation and energy spectrum shift, are essentially captured by the 
suggested relativistic extension of our quantum theory. 

7. Conclusions 

In this paper, a new mathematical formulation of quantum mechanics has 
been suggested in terms of signed Newtonian field-less particles which can 
be seen as a generalization and physical interpretation of the time-dependent 
Wigner Monte Carlo method [H] extended to a non-discretized phase-space and 
a non-finite domain. The suggested formulation consists of a set of three pos¬ 
tulates which completely describe the time evolution of quantum systems. The 
predictions made by this theory are exactly the same as those made by other 
standard approaches. In fact, we have shown that simulating a quantum system 
by repeatedly applying the operator S is perfectly equivalent to solve the time- 
dependent Wigner equation and, consequently, the time-dependent Schrodinger 
equation. In order to show the validity and applicability of the suggested for¬ 
mulation, we have simulated different systems where quantum tunnelling effects 
are prominent. Finally, we have proposed a relativistic extension of the second 
postulate in terms of signed particles moving on space-time geodesics and we 
have shown how typical relativistic effects, such as time dilation and shift of 
energy spectrum, are predicted correctly. This particular example underlines 
how the new formulation provides a novel significantly easy-to-extend quantum 
theory. 

To conlude, one should note that many aspects of this peculiar and new 
formulation still remain to be explored. For instance, it would be interesting 
to study how the introduction of a finite time step At affects the simulation of 
a quantum system (as it may strongly depends on the presence of an external 
potential). Likewise, the frequency of annihilation steps to apply in a simu¬ 
lation remains an open problem which needs to be investigated. It would be 
also important to study the mathematical properties (and in particular the con¬ 
vergence) of the newly defined momentum integral Q. Moreover, it would be 
useful to investigate the relation between this new formulation and others such 
as the ones proposed by D. Bohm and R. Feynman (as many similarities can 
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be found already). It is also clear that the classical limit presented in this work 
may not be entirely satisfactory from a mathematical/theoretical perspective as 
it introduces several strong assumptions which may need to be investigated in 
further details. Another interesting point is to apply this method to bigger and 
more realistic quantum systems occurring in real life applications and see how it 
performs compared to other more standard approaches. This will be the topic 
of next-future papers. Despite all these pending questions, the signed particle 
formulation remains a new instrument to explore quantum mechanics, offering 
a unique and alternative new point of view. 

8. Computational aspects 

The simulator used to obtain the results presented in this paper is a modified 
version of Archimedes, the GNU package for the simulation of carrier transport 
in semiconductor devices Hi. This code was first released in 2005, and, since 
then, users have been able to download the source code under the GNU Public 
License (GPL). Many features have been introduced in this package. In this 
particular project, our aim has been to develop a full quantum time-dependent 
nanodevice simulator including phonon scattering effects. The code is entirely 
developed in C and optimized to get the best performance from the hardware. 
It can run on parallel machines using the OpenMP standard library. The results 
of the present version are posted on the nano-archimedes website, dedicated to 
the simulation of quantum systems. |19j . 

The results have been obtained using the HPG cluster deployed at the Insti¬ 
tute of Information and Gommunication Technologies of the Bulgarian Academy 
of Sciences. This cluster consists of two racks which contain HP Gluster Plat¬ 
form Express 7000 enclosures with 36 blades BL 280c with dual Intel Xeon 
X5560 @ 2.8 Ghz (total 576 cores), 24 GB RAM per blade. There are 8 storage 
and management controlling nodes 8 HP DL 380 G6 with dual Intel X5560 @ 
2.8 Ghz and 32 GB RAM. All these servers are interconnected via non-blocking 
DDR Infiniband interconnect at 20Gbps line speed. The theoretical peak per¬ 
formance is 3.23 Tflops. 
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Figure 1: Quantum tunnelling experiment in phase-space: an initially Gaussian wave-packet 
is evolving in free space. The wave packet is represented by a set of signed particles (initially 
all positive) which is evolved in time by repeatedly applying the operator S. This plot shows 
the simulation at times (from left to right, from top to bottom) 0 fs, 10 fs, 20 fs, 30 fs, 40 fs, 
50 fs, 60 fs and 70fs respectively. The (red) liji 0 represents the middle of the spatial domain. 
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Figure 2: Quantum tunnelling experiment in phase-space: an initially Gaussian wave-packet is 
directed towards a potential barrier equal to 0.10 eV. The wave packet is represented by a set 
of signed particles (initially all positive) which is evolved in time by repeatedly applying the 
operator S. This plot shows the simulation at times (from left to right, from top to bottom) 
0 fs, 10 fs, 20 fs, 30 fs, 40 fs, 50 fs, 60 fs ar^TOfs respectively. One notes the appearance 
of negative values (dark blue) during the evolution. The position of the potential barrier is 
schematically represented by the (red) line in the middle of the domain. 


















































Figure 3: Quantum tunnelling experiment in phase-space: an initially Gaussian wave-packet is 
directed towards a potential barrier equal to 0.30 eV. The wave packet is represented by a set 
of signed particles (initially all positive) which is evolved in time by repeatedly applying the 
operator S. This plot shows the simulation at times (from left to right, from top to bottom) 
0 fs, 10 fs, 20 fs, 30 fs, 40 fs, 50 fs, 60 fs ar^TOfs respectively. One notes the appearance 
of negative values (dark blue) during the evolution. The position of the potential barrier is 
schematically represented by the (red) line in the middle of the domain. 
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Figure 4: Proper time of an initially Gaussian wave-packet in the non-relativistic case, dashed 
(blue) curve, and relativistic case, continuous (red) curve, respectively. The phenomenon of 
time dilation, typical of relativistic theories, is clearly observable. 
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Figure 5: Induced dipole of an initially Gaussian wave-packet in the non-relativistic case, 
dashed (blue) curve, and relativistic case, continuous (red) curve, respectively. A shift in the 
energy spectrum, typical of relativistic systems, is clearly visible. 
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